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Abstract 

In this paper we discuss results relevant to 3D Double-Side Double Type Column (3D-DDTC) pixel sensors fabricated at FBK 
(Trento, Italy) and oriented to the ATLAS upgrade. Some assemblies of these sensors featuring different columnar electrode 
configurations (2, 3, or 4 columns per pixel) and coupled to the ATLAS FEB read-out chip were irradiated up to large proton 
fluences and tested in laboratory with radioactive sources. In spite of the non optimized columnar electrode overlap, sensors exhibit 
reasonably good charge collection properties up to an irradiation fluence of 2x 10^^ neqCm"^, while requiring bias voltages in the 
order of 100 V. Sensor operation is further investigated by means of TCAD simulations which can eff'ectively explain the basic 
' mechanisms responsible for charge loss after irradiation. 
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1. Introduction 

The fast increase in luminosity in the modern High Energy 
Physics (HEP) experiments is pushing the research in the field 
of silicon radiation detectors to new challenging frontiers. Due 
to the high radiation doses foreseen for the inner tracking lay- 
ers, radiation hard detectors must be designed and tested in or- 
der to provide reliable particle detection up to fluences in the 
order of 10^^ 1-MeV equivalent neutrons per square centimeter 
(neqcm"^). At the same time these devices must be fast in term 
of charge collection time and less power consuming than the 
older ones. For these reasons several R&D projects in the field 
of silicon radiation detectors have been launched in the past 
years, mostly focusing on the upgrades of the experiments at 
the Large Hadron Collider (LHC) at CERN, Geneva, Switzer- 
land fri. 

The main macroscopic consequences of radiation-induced de- 
fects in the detector bulk are: (i) changes in eff'ective dop- 
ing concentration, mainly with introduction of acceptor-like 
defects that lead to an increase of the full depletion voltage, 
(ii) higher leakage currents due to the creation of genera- 
tion/recombination centers, and (iii) decrease of the charge col- 
lection efficiency due to carrier trapping O. The overall con- 
sequence of this damage is a strong reduction in the signal to 
noise ratio that can severely reduce the tracking capabilities. To 
counteract these eff'ects diff'erent strategies are possible |3|: i) 
material engineering, i.e., using as a substrate either non stan- 
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dard silicon (e.g.. Magnetic Czochralski, epitaxial, etc.) or di- 
amond, which are intrinsically more resistant to radiation dam- 
age; ii) device engineering, which consists in designing detec- 
tors with geometrical configurations that allow for lower signal 
degradation after irradiation. One of the most promising ap- 
proaches to achieve radiation hard silicon detectors is the so- 
called 3D-architecture proposed by Parker and collaborators in 
the mid '90s |4|. In 3D detectors the electrodes have a colum- 
nar shape and are etched perpendicularly to the wafer surface, 
penetrating the entire sensor thickness. In this way the distance 
between electrodes is not bound to the thickness of the wafer 
(which is the case for standard planar silicon detectors) but can 
be optimized to suit performance requirements. Thanks to this 
characteristic the distance between the electrodes is decoupled 
from the active volumes thickness. As a consequence, low oper- 
ating voltages (less than 10 V before irradiation, at most 200 V 
after irradiation), fast response times, and strong reduction of 
charge trapping eff'ects after irradiation are obtained Q. An- 
other important feature deriving from 3D detector technology is 
the active edge, which consists of a trench electrode termination 
allowing for a good sensitivity up to a few microns aways from 
the physical edge of the sensors. As a result a more efficient 
area coverege on wide surfaces and lower material budget are 
obtained 01 . While active edge is an intrinsic option for 3D de- 
tectors, it can also be implemented in planar sensors, although 
with a major process complication |7 1. Besides all these advan- 
tages, 3D detectors have some disadvantages: in particular, the 
fabrication process is more complicated than a standard silicon 
detector process, the capacitance is higher and their response is 
not completely uniform because of the electrodes, that are not 
fully efficient, and of the presence of some low field regions. In 
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order to develop 3D silicon detectors for the ATLAS upgrade 
the so-called ATLAS 3D sensor collaboration 1 8 1 was formed, 
involving many research centers and institutes from all over the 
world. Among the technological approaches considered for 3D 
fabrication, besides the original one developed at Stanford |9|, 
there are also simplified architectural implementations. One of 
them, relevant to the detectors considered in this work, is the so- 
called 3D Double sided Double Type Column (3D-DDTC) con- 
cept, independently proposed by FBK, Trento, Italy 1 10] and by 
CNM, Barcelona, Spain |11| with the aim of reducing process 
complexity in view of medium volume productions. One of 
the main advantages of this approach is that it does not use a 
support wafer, thus avoiding the related steps of wafer bond- 
ing and final wafer removal. Moreover, in 3D-DDTC detectors 
the substrate bias can be applied from the back side, making 
these sensors compatible with standard planar sensors and eas- 
ing the detector assembly within a tracking system. Columns 
are etched from both wafer sides (n^ from the top, p^ from the 
bottom) and do not pass through the entire wafer thickness, so 
they only partially overlap. From TCAD simulations ifTOl , it 
was predicted that the performance of 3D-DDTC detectors is 
comparable to that of standard 3D detectors if the column over- 
lap is a significant fraction of the wafer thickness, whereas it 
can be degraded if column thicknesses are not optimized, which 
is the case for the first prototypes fabricated at FBK and consid- 
ered in this paper. Therefore, the radiation hardness should be 
carefully studied in order to obtain useful information for the 
design and technology optimization. The FBK devices were 
previously tested both in laboratory 1 12] and in beam tests at 
CERN in pre-irradiation conditions 1 13], obtaining very good 
results. In order to study their radiation hardness, diff'erent ir- 
radiation campaigns were conducted, and irradiated detectors 
were measured again in a test beam at CERN 1 14] and in labo- 
ratory. 

In this paper we report on selected results from functional char- 
acterization with radioactive sources conducted in laboratory 
on these 3D-DDTC detectors. Numerical simulations are also 
used to gain better insight into experimental results. The paper 
is organized as follows: Section 2 gives a brief description of 
devices under test and summarizes the two proton irradiation 
campaigns; Section 3 describes and discusses post-irradiation 
measurement results also comparing them with pre-irradiation 
results; Section 4 reports numerical simulation results and com- 
pares them to the measurements. Conclusions follow. 

2. Experimental 

2. 1 . Sensor description 

The sensors under test are 3D-DDTC detectors fabricated at 
FBK on 4", 200 yum thick, p-type, FZ silicon wafers. As pre- 
viously stated, in these devices columns are etched from oppo- 
site sides of the wafer and are not completely passing through 
the silicon bulk (91. At the time of fabrication of these sen- 
sors (2008) the Deep Reactive Ion Etching (DRIE) equipment 
was not yet available at FBK so the etching of the holes was 
commissioned to an external company (IBS, France) and prob- 
lems related to the calibration of this step led to an un-even 



column depth (see Fig. [T] (left)). This problem translated into a 
relatively small column overlap, in the order of 90 yum, which 
of course is not ideal and could aff'ect the device performance, 
especially after irradiation. The nominal column diameter is 
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Figure 1 : (left) Schematic cross-section of the sensors and (right) different pixel 
configurations. 

lOyum. The surface insulation between n^ electrodes on the 
front side is achieved by combined p-spray and p-stop implants 
1 15], whereas all ohmic columns are shorted on the back side 
by uniform p^ diffiision and metal. 

Devices under test are pixel detectors compatible with the FEI3 
ATLAS read-out chip 1 16|, and feature various layout options 
diff'ering in the number of columns per pixel: 2E-type (two 
n^ columns per pixel), 3E-type (three n^ columns per pixel) 
and 4E-type (four n^ columns per pixel). The number of 
ohmic columns per pixel also changes accordingly, and the 
inter-electrode distances are 103 yum, 71.2yum, and 56yum, re- 
spectively (see Fig. [T] (right)). Detectors were bump bonded to 
FEB ATLAS read-out chips at SELEX (13, Rome, Italy. The 
FEB chip was designed with radiation tolerant layout rules in a 
0.25 yum CMOS technology. The maximum radiation dose that 
the chip can withstand is in the order of 50Mrad. Detectors 
were designed to exactly match the geometry of the readout 
channels of the chip. The front-end chip presents 2880 chan- 
nels which are arranged in a matrix of 160 rows per 18 columns, 
with a pixel size of 400 yum x 50 yum. Each channel is composed 
by an analog and a digital part. The analog part integrates the 
sensor output current by means of a charge sensitive preampli- 
fier with constant current discharge, thus yielding a triangular 
pulse shape, which is fed to a discriminator along with a pre- 
set threshold. As a result, the width of the discriminator output 
signal, i.e., the Time Over Threshold (TOT) expressed in units 
of 40 MHz clock, is in first approximation proportional to the 
collected charge. A detailed explanation of the FEB operation 
can be found in |[T6ll . 

2.2. Proton irradiation campaigns 

In order to study the radiation tolerance of these sensors, five 
devices were irradiated to two diff'erent proton fluences at two 
facilities: three samples at the Karlsruhe Institute of Technol- 
ogy (KIT) with a 25 MeV proton beam up to 1 x lO^^neqCm"^ 
and two samples at CERN PS with a 24GeV/c proton beam 
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Module ID 


Sensor Type 


Fluence [neqCm ^] 


Particle 


Facility 


A 


2E 


1x10^^ 


25 MeV proton 


KIT 


B 


3E 


1x10^5 


25 MeV proton 


KIT 


C 


4E 


1x10^5 


25 MeV proton 


KIT 


D 


2E 


2x10^5 


24 GeV proton 


CERN PS 


E 


4E 


2x10^5 


24 GeV proton 


CERN PS 



Table 1: Overview of irradiated assemblies. 



up to 2x lO^^neqCm"^ (see Table[T]). The proton fluences were 
scaled to IMeV equivalent neutrons per square centimeter 
(neqCm"^) using the NIEL hypothesis with hardness factors of 
1.85 and 0.62 for 25 MeV and 24 GeV protons, respectively. 
The uncertainty in the irradiation fluences is lower than 10%. 
After the irradiation, the detectors were cooled to prevent an- 
nealing. However, the detectors had to be kept at room temper- 
ature for a short time during handling and measurement setup, 
so that they experienced some annealing (which was of frac- 
tions of an hour). 

3. Functional lab-test measurements 

Results from the electrical and functional characterization 
before irradiation of pixel detectors here considered are re- 
ported in 1 12], |18|. In this paper we will focus on mea- 
surements performed on irradiated devices and results of non- 
irradiated detectors will be recalled for a direct comparison 
(see a summary in Table [2]). The sample characterization has 
been carried out by measuring: leakage current versus bias 
voltages, threshold, noise and response to radioactive y- and 
yS-sources. All measurements on irradiated devices were per- 
formed at -20 °C in order to reduce leakage current and avoid 
reverse annealing efl'ects. 
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Figure 2: Overview of I-V curves from irradiated assemblies. 

Before irradiation detectors had leakage currents of a few hun- 
dred nA at room temperature and breakdown voltages between 
60 V and 70 V. After irradiation, the leakage currents and the 



breakdown voltages shift to higher values due to the displace- 
ment damage of the silicon and to the increase of the oxide 
charge concentration caused by ionizing radiation. It was mea- 
sured that the leakage current increased by several orders of 
magnitude, whereas, the breakdown voltages were generally 
increased, sometimes well above lOOV, although results are 
not completely uniform. Breakdown voltage improvement af- 
ter irradiation is typical of p-spray isolated planar structures 
ifTSl . These 3D-DDTC sensors confirm the expectations. Fig.|2] 
shows an overview of leakage current versus bias voltage: the 
current values and trends are in agreement with expectations, 
with diff'erent current levels corresponding to diff'erent irradi- 
ation fluences. Assemblies A, B and C were irradiated at 
1x10^^ neqCm"^ and show similar current values after irradi- 
ation (reaching about 25 jiA). A clear saturation at about 100 V 
can be observed for samples A and B, whereas sample C suf- 
fers from an earlier breakdown at about 65 V. Assemblies D 
and E were irradiated at 2x 10^^ neqCm"^ and reach currents of 
35 to 40yuA (not yet fully saturated) before breakdown. The 
agreement between the leakage current damage constant (a ^ 
5x10"^^ A/cm) extracted from these measurements and the 
generally accepted value of 4 x 10"^^ A/cm |2| is good enough 
considering the uncertainties in the irradiation fluence, anneal- 
ing conditions and temperature. 

In order to evaluate the performance of the assemblies under 
test (sensor and front-end electronics), the system was cali- 
brated aiming at a ToT of 60 units for a charge of 20 ke and 
a threshold of 3.2 ke. Threshold and noise measurements have 
been performed on each pixel based on S -curve fit function 
|[T9l . Table [3] summarizes the values of threshold and noise for 
the irradiated assemblies. Compared to values measured before 
irradiation (cf. Table |2]), both threshold average and dispersion 
are very similar. Also noise values, expressed as equivalent 
noise charge (ENC), are only slightly diff'erent than before ir- 
radiation. As an example. Fig. [3] shows the ENC as a function 
of the voltage for sensors A, B, and C. The ENC curves are al- 
most flat because, after irradiation to such a large fluence, the 
capacitance, which is the main factor for noise, is almost con- 
stant with bias due to the very high resistivity of the substrate 
II2TII . Difl'erent noise levels are observed, in agreement with the 
pre-irradiation case, because of the difl'erent capacitances char- 
acterizing the difl'erent column configurations (cf. Table [2]). A 
direct comparison in the noise values before and after irradia- 
tion is difficult because of the different temperatures (-h20 °C 
vs -20 °C): on one hand, the increase of leakage current in ir- 
radiated sensors would justify an increase in the noise in the 
order of tens of electrons rms 1201 , on the other hand this in- 
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Sensor type 


^leakage 


[nA] 


^breakdown [V] 


Threshold [e] 


ENC [e] 


Am 60-keV peak [ke] 


^"Sr MPV peak [ke] 


2E 


250- 


300 


60-70 


3240 ± 54 


197 + 9 


14.71 


15.35 


3E 


250- 


300 


60-70 


3274 + 51 


206 + 8 


14.50 


15.25 


4E 


250- 


300 


60-70 


3297 ± 56 


227 + 8 


14.37 


15.25 



Table 2: Summary of typical results obtained from non-irradiated sensors. Except for Vbreakdown. the parameters are obtained at a bias voltage of 35 V. 



Module ID 


Sensor type 


Threshold [e] 


ENC [e] 


Vtest [V] 


A 


2E 


3261+73 


204.1 + 11.6 


120 


B 


3E 


3158 + 140 


233.2+ 13.4 


100 


C 


4E 


3267 + 78 


232.4+12.7 


60 


D 


2E 


2950 + 206 


119.4 + 31.5 


125 


E 


4E 


3307 + 115 


189.6+19.0 


80 



Table 3: Threshold and noise values measured from irradiated assemblies. 
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Figure 3: Equivalent noise charge as a function of bias voltages for irradiated 
assemblies A, B, and C. 



crease might well be compensated by the much lower tempera- 
ture used for measurements after irradiation. Limiting the com- 
parison to the irradiated samples, it might be puzzling that those 
irradiated at 2x 10^^ neqCm"^ (samples D and E) show lower 
noise than those irradiated at 1 x 10^^ ngqCm"^ (samples A, B, 
and C), but it should be noted that this difference might be due 
to the radiation damage to the electronics rather than to the sen- 
sors: in fact, for the 24-GeV proton irradiation a Total Ionizing 
Dose (TID) of about 94 Mrad has been estimated, whereas the 
estimated TID is 144 Mrad for 25-MeV proton irradiation. Con- 
sidering that the letter TID is about 35% higher than the former 
and about three times the one for which the FEB chip has been 
designed (50 Mrad), it could be responsible for a degradation in 
the chip noise performance. 

Charge collection mechanisms were studied by means of an 
^"^^ Am y-source and a ^^Sr yS-source, comparing the results with 
those obtained before irradiation (cf. Table [2]). For ^"^^ Am mea- 
surements, the self-triggering capabilities of the system were 



exploited, whereas for ^^Sr measurements the trigger was taken 
from a scintillator placed behind the devices under test. 
Tableg summarizes the results of ^"^^Am y-source measure- 
ments for the three sensors irradiated at 1 x 10^^ neqCm"^ (data 
are not available for samples D and E). Before irradiation (cf. 
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Figure 4: Mean peak (60-keV) of the collected charge for Am source mea- 
surements as a function of reverse bias voltage in all irradiated samples. 

Table [2]) the mean of the fitted Gaussian of about 14.5 ke with- 
out clustering had been obtained with an applied reverse bias of 
35 V. This charge corresponds to the energy peak at 60 keV, in 
agreement with theoretical expectations within the uncertainty 
due to the calibration process, which was estimated to be in 
the order of 10-15 %. This uncertainty is due to the Time over 
Threshold (ToT) to charge calibration process which contains 
several components and it estimated to be about mentioned per- 
centage. In a first approximation, the 10-15 % of uncertainty 
is related to the uncertainty of the capacitance used to tune the 
first stage of the EE in the calibration process. Detailed expla- 
nation can be found in 1 16 |. 

Looking at Table[4j it is possible to notice that after irradia- 
tion, samples B and C yield a value of collected charge com- 
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Module ID 


Sensor type 


Am 60-keV peak [ke] 


Vtest [V] 


A 


2E 


12.8 


80 


B 


3E 


14.5 


100 


C 


4E 


14.3 


65 



Table 4: Collected charge from the assemblies irradiated at 1 x 10^^ neqcm ^ exposed to Am source. 



Module ID 


Sensor type 


^^Sr MPV peak [ke] 


Vtest [V] 


A 


2E 


10.2 


120 


B 


3E 


11.9 


100 


C 


4E 


10.9 


60 


D 


4E 


9.0 


125 


E 


4E 


9.8 


80 



Table 5: Collected charge from all irradiated assemblies exposed to ^^Sr source. 



parable to the value before irradiation, whereas sample A (2E- 
type) shows lower collected charge. This can be ascribed to the 
longer distance between the n^ and p^ electrodes in sample A, 
which results in higher trapping probability. The importance of 
the column geometry is confirmed by results in Fig.|4] which 
shows the mean charge values as a function of the bias voltage 
for the three samples. Sensor C (4E-type) reaches roughly the 
same maximum charge as sensor B (3E-type), but at a lower 
voltage, owing to a lower distance between electrodes. Even 
higher voltage is necessary to have a significant charge collec- 
tion in sensor A, which also tends to saturate at a lower charge 
value. An example of ^^^Am spectrum measured with sample C 
is shown in Fig. [5] The main peak corresponding to the 60-keV 
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Figure 5: Am spectrum measured with a 4E-type detector (sample C) irradi- 
ated with protons at 1 x 10^^ neqcm"^ and reverse biased at 65 V. 

photons is distorted with a tail towards lower charge values, 
that can be attributed to the absorption of photons in regions 
with lower electric field, for which the charge collection is not 
fully efficient. 

Table [5] summarizes yS-source measurements for all irradiated 
detectors biased at the optimal voltage. The devices have been 
measured at diff'erent bias voltages, always stopping at a volt- 



age lower than the breakdown voltage. Fig. 6(a) shows the pulse 
height spectrum in response to a ^^Sr yS- source in the 4E sensor 
irradiated with protons at 1 x 10^^ neqCm"^ (sample C). The dis- 
tribution is related to all clusters, and charge values in excess 
of 7 ke have been fitted with a Landau function, which is also 
shown in the figure. The most probable value (MPV) of the col- 
lected charge is 10.65 ke, whereas it was about 15 ke before irra- 
diation (cf. Table[2]). Thus, the degradation in charge collection 
is more pronounced than that observed for the Am y-source 
measurements. This is generally the case for all the irradiated 
assemblies, as can be seen from the values in Table (5] The rea- 
son for this behavior is related to the diff'erent ways charge is 
generated inside the sensors: while every y-ray in silicon can 
be absorbed at diff'erent depths and the released charge is lo- 
calized in a region of a few cubic micrometers, /3 particles tra- 
verse the whole bulk thickness releasing charge along the entire 
track. This means that, if the device is not fully depleted, part 
of the charge released by a particle will be generated in a non 
depleted region and therefore almost completely lost (due to re- 
combination). On the contrary, for all the y-rays absorbed in the 
region where columns overlap the generated charge will be col- 
lected also after irradiation, except for the case where the inter- 
electrode distance is too large, because of trapping. Assemblies 
irradiated at 1 x 10^^ ngqCm"^ perform reasonably well, consid- 
ering that the applied voltage is high enough to achieve lateral 
depletion between the columns. However, it should be stressed 
that a non negligible fraction of the active volume at the bot- 
tom is not depleted, thus causing charge loss (this eff'ect will 
be better explained in the following with the aid of TCAD sim- 
ulations). Samples irradiated at 2x 10^^ neqCm"^ show further 
degradation of the collected charge with respect to those irradi- 
ated at 1 X 10^^ neqcm"^, due to stronger trapping eff'ects. Nev- 
ertheless, the collected charge is still reasonably good because 
of the higher bias voltage that could be applied to these sensors 
before they reached breakdown, which allowed lateral deple- 
tion between columns to be achieved both for the 2E and the 
4E configurations. Again, the 4E sample collects more charge 
than the 2E one at a lower voltage owing to the shorter distance 
between the electrodes. The cluster size distribution for sample 
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Figure 6: Results for tests with ^*^Sr source on a 4E-type detector (sample C) 
irradiated with protons at 1 x 10^^ neqcm"^ and biased at 60 V: (a) pulse height 
distribution, and|(b)|cluster size distribution. 



C is also shown in Fig. |6(b)| Most events (-80 %) are cluster 
size 1, the remaining part (~20 %) are cluster size 2. This dis- 
tribution can be explained taking into account that ft particles 
reaching the sensor are not all perpendicular to the surface due 
to the geometry of the collimator. The angle distribution was 
found to span from 90 ° to 84.29 °. Although in 3D sensors the 
electric field distribution provides a sort of self-shielding eff'ect 
in each pixel, in case particle hits are not perpendicular to the 
surface the probability of charge sharing between two adjacent 
pixels becomes non negligible. 

The performances of some of the irradiated devices considered 
in this paper were also tested with a 120 GeV/c tt^ beam at 
CERN SPS in June 2010. The aim of this test was to study 
tracking efficiency, charge sharing and cluster size. A detailed 
description of the obtained results can be found in |14|. It 
should be noticed that both the values of the collected charge 
and the cluster size distributions are in very good agreement 
with those measured in laboratory. 

The MPV of the collected charge reported in Table |5] are rele- 
vant to the highest possible bias voltage. The variation of the 



MPV with the bias voltage is also of interest to better under- 
stand the behavior of these sensors. Fig. [7] summarizes data 
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Figure 7: Most Probable Value (MPV) of the collected charge for ^^Sr source 
tests as a function of reverse bias voltage in all irradiated samples. 

for all the irradiated assemblies. As expected, the collected 
charge increases with the applied voltage in agreement with 
the larger depleted volume within the sensors. It should also 
mentioned that a loss of efl&ciency associated to the increasing 
number of electrodes was measured in the CERN SPS beam- 
tests as is reported in |[T4ll . Apart from sensor A, for which 
trapping effects are more severe as also observed with y-rays, 
charge values are not yet completely saturated in the consid- 
ered voltage range. As expected, at the same bias voltage the 
collected charge increases with the number of electrodes. The 
same value of collected charge is observed at a lower voltage 
in 4E samples with respect to 3E and 2E ones. Samples irra- 
diated at 2x 10^^ neqcm"^ of course need a higher bias voltage 
to reach a sufficient charge collection efficiency, and maximum 
values are slightly lower due to trapping. 



4. TCAD simulations 

In order to better understand the experimental results, numer- 
ical simulations were performed with TCAD tools from Synop- 
sys |22|. Simulations were focused on devices measured with 
^^Sr yS-source and irradiated at at 1 x 10^^ ngqCm"^, in order to 
limit the computational load and reduce the simulation time. 
Simulations were performed at different bias voltages. More- 
over, the two extreme values of particle incidence angle were 
tested (90° and 84.29 °). To this purpose, since the probability 
of charge sharing is strictly related to the inclination of the par- 
ticle, two different structures were simulated see Fig. 8, which 
refers to the case of the 4E devices): (i) a single cell including 
one n^ column and one p^ column, for a 90 ° hit angle (Fig. 8(a) 



and (ii) two adjacent cells including two n^ electrodes (related 
to adjacent pixels) and one p^ column, for a 84.29 ° hit angle 



(Fig. 8(b) ). Moreover, to have a better understanding of the be- 
havior of the detectors, different particle hit points were tested 
and the results were combined as explained below. 
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Figure 8: Front view of 4E simulated structures: |(a)|st ructure for cluster size 1 
simulations with different particle hit points, and |(b)| structure for cluster size 
2 simulations with different particle hit and arrival points (track angle 84.29 °). 
Similar simulation domains choice apply to the 2E and 3E devices. 



A particle hit was simulated using the Heavy-Ion model avail- 
able in the simulator: the released charge was 80 electron-hole 
pairs per micron and the spatial distribution was Gaussian in a 
region of one micrometer diameter around the track. The radi- 
ation damage was modeled using the Perugia trap model |23 1 
modified as described in L24J . 

Transient simulations were performed from to 100 ns; the 
high leakage current was subtracted from the output current and 
a numerical integration was performed in order to extract the to- 
tal collected charge for each simulation. In order to decide if an 
event was good or not, the collected charge after 20 ns was ob- 
served and, if the value was above the threshold (3200 e), the hit 
was considered valid and the value at which the integral satu- 
rated was taken as the total charge collected in that event. In the 
charge sharing simulations, the same procedure was adopted for 
both n^ columns. 

Since measurements showed diff'erent distributions for cluster 
size 1 and cluster size 2 events, simulations were combined us- 
ing the following equation: 



Figure 9: Comparison between measured and simulated collected charge as a 
function of reverse bias voltage for sensors at 1 x 10^^ neqcm"^. 



where Qcsi is the total collected charge for cluster size 1 
events, Qcsi is the total collected charge for cluster size 2 
events (the sum of the charge collected on two adjacent pixels) 
and CSl% and CS2% are the percentages of cluster size 1 
and cluster size 2 events obtained from the measurements. In 
particular the 2E and 3E devices showed 90% of cluster size 1 
events and 10% of cluster size 2 events while these percentages 
became 80% and 20% for the 4E device. Simulation results 
were compared with measured Most Probable Values and a 
very good agreement was found (see Fig.|9]), thus confirming 
that the most important physical mechanisms are properly 
modeled in the simulations. In order to gain a better insight into 
the charge collection process, several cuts were extracted from 
the simulations results showing diff'erent electrical quantities, 
among them the electric field, Shockley-Read-Hall (SRH) 
recombination, electron density, and hole density for diff'erent 
bias voltages and diff'erent depths inside the device. Let's 
consider as an example the 4E device. Two voltages (20 V and 
60 V) were chosen as representative of low-bias and high-bias 
conditions, and the cuts were extracted from two separate re- 
gions, one in the part of the sensor where columnar electrodes 
overlap (z = 50yum) and the other in the non-overlap region 
(z = 150yum). This was done as these two regions of the device 
behave completely diff'erently, especially after irradiation. 
The above-mentioned quantities, relevant to a time instant 
0.5 ns after the particle hit the device perpendicularly to the 



wafer surface (point E in Fig. 8(a) ), are shown in Figs. 10 - 12 



Looking at Fig.[TO| it is clear that for both bias voltages the 
electric field is strong only in the overlap region (with a much 
higher absolute value for 60 V bias), whereas it is very weak in 
the bottom part of the sensor. 



From Fig. 1 1 it is possible to confirm that SRH recombination 
is negligible in the overlap region, whereas it significantly 
aff'ects those carriers generated by the impinging particle in the 
non-overlap region. 



Qall-cluster,sim - CS 1% X QcSl + 2% X QcS2 (1) 
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Figure 10: 2d cuts in the 4E sensor showing the electric field at two diff'erent depths (Z) and two bias voltages (Vb): (A) Z = 50//m, Vb = 20V; (B) Z = 50//m, 
Vb = 60 V; (C) Z = 150 /um, Vb = 20 V; (D) Z = 150 /um, Vb = 60 V 
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Figure 1 1 : 2d cuts in the 4E sensor showing 

Vb = 60 V; (C) Z = 150 jim, Vb = 20 V; (D) Z = 150 //m, Vb = 60 V 
surface (point E in|8(a)). 



showing the SRH recombination at two diff'erent depths (Z) and two bias voltages (Vb): (A) Z = 50 /urn, Vb = 20 V; (B) Z = 50 /urn, 
^ \i- ^T^\ v - i ,,rr, \l - /^n \i D^ta are relevant to a time instant 0.5 ns after a particle hit the device perpendicularly to the wafer 
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Figure 12: 2d cuts in the 4E sensor showing the electron density at two different depths (Z) and two bias voltages (Vb): (A) Z = 50 jum, Vb = 20 V; (B) Z = 50 jum, 
Vb = 60 V; (C) Z = 150 /urn, Vb = 20 V; (D) Z = 150 /urn, Vb = 60 V Data are relevant to a time instant 0.5 ns after a particle hit the device perpendicularly to the wafer 
surface (point E in|8(a)). 
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Figure 13: Simulated output currents in the 4E sensor at 20 V and 60 V bias |(a)| and corresponding time integrals |(b)] 
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Both the distributions of the electric field distribution and the 
SRH recombination can be explained considering the way the 
depletion region extends inside this type of devices: full deple- 
tion occurs at relatively low bias voltages in the overlap region, 
whereas much larger voltages would be required to fully de- 
plete the bottom part of the device. 

Given that part of the charge generated by the impinging par- 
ticle is lost because of recombination, it is important to under- 
stand how the rest of the carriers move inside the device and 



how they contribute to the signal formation. Fig. 12 shows the 
electrons density distributions for the biasing conditions, depths 
and time instant. From the electron density point of view, it is 
clear that carriers generated in the top part of the detector are 
collected almost immediately owing to the strong electric field. 
The eff'ect of the bias voltage is clearly visible because the elec- 
trons are collected faster at 60 V than at 20 V in this region. The 
situation in the bottom part of the device is diff'erent: the elec- 
tric field is much weaker and carriers must drift (or even diff'use 
in non depleted regions) both laterally and vertically toward the 
tip of the n^ column. Since distances are longer and the field is 
weaker, carriers in this region will require more time to be col- 
lected, drastically increasing the probability of being trapped 
or recombining. The minor changes in the bottom part of the 
device between 20 V and 60 V indicate that the device has not 
reached an optimal working point. The previous considerations 
can also be applied to what happens for the hole density (not 
shown). Again, holes are collected much faster in the overlap 
region than in the bottom region. 

The above analysis is in good agreement with simulation re- 



sults in terms of transient output currents (Fig. 13(a)) and their 



time integrals (Fig. | 13(b)] ) for the two considered biases (20 V 
and 60 V). As already mentioned, for higher bias voltages, elec- 
trons and holes in the overlap region suff'er from less trapping 
and this translates into a pulse with higher amplitude and faster 
evolution thanks to the higher electric field (both electrons and 
holes generated in the overlap region concur in the formation 
of the fast peak). Apart from the diff'erence in the peak ampli- 
tude and time, a change of the shape of the pulse tails can be 
observed: at 60 V of bias, the pulse tail extends to a slightly 
longer time, a fact that is also confirmed by the delayed satura- 
tion of the corresponding time integral curve in Fig. | 13(b)] The 
change of the shape of the pulse tails is related to the fact that 
at 60 V of bias, the bottom part of the detector suff'ers from less 
recombination and the charge generated there also contributes 
to the output signal. Since the field in this region is weaker 
than in the top, the drift of the carriers will be slower, hence the 
presence of the long tail. Similar considerations apply to the 
other types of sensors, that are not reported here for the sake of 
conciseness. 

5. Conclusion 

We have reported on selected results from the functional 
characterization of irradiated 3D-DDTC pixel sensors fabri- 
cated at FBK. Sensors with diff'erent pixel configurations have 
been assembled with the ATLAS FEB read-out chip and ir- 
radiated with protons up to very large fluences. Experimen- 



tal results from measurements carried out in the laboratory 
have been discussed and compared to those obtained before 
irradiation. As expected, the breakdown voltage was gener- 
ally increased by a few tens of Volts with respect to the pre- 
irradiation values, allowing to operate the sensors at a bias volt- 
age high enough to achieve at least lateral depletion between 
the columns. The peak of the collected charge in response 
to y-rays from an ^"^^Am source is not significantly degraded 
with respect to the values obtained before irradiation, provided 
that the bias voltage is high enough to eff'ectively counteract 
trapping in the top region of the devices where columns over- 
lap. On the contrary, the most probable value of the collected 
charge for tests with yS-particles from a ^^Sr source is found to 
be more sensitive to trapping, due to the lower charge collec- 
tion efficiency from the bottom region of the devices, where 
columns do not overlap, in good agreement with TCAD simu- 
lations. Nevertheless, the signal efficiency, defined as the ratio 
of the MPV of the collected charge after irradiation and before 
irradiation, remains at acceptable levels: for the best sample ir- 
radiated at 1 X 10^^ neqCm"^ it is about 76 % at 100 V, and for 
the best sample irradiated at 2 x 10^^ neqCm"^ it is about 64 % at 
120 V. These results are very encouraging since the tested sam- 
ples have a rather short column overlap, less than one half of 
the sensor thickness. Hence, there is still wide room for perfor- 
mance improvement by etching deeper junction column depths. 
New 3D-DDTC sensors fabricated at FBK have indeed passing- 
through columns and are expected to be very radiation hard. 
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